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ABSTRACT 

We study the magnetosphere of a slowly rotating magnetized neutron star subject to 
toroidal oscillations in the relativistic regime. Under the assumption of a zero inclination 
angle between the magnetic moment and the angular momentum of the star, we analyze the 
Goldreich- Julian charge density and derive a second-order differential equation for the elec- 
trostatic potential. The analytical solution of this equation in the polar cap region of the mag- 
netosphere shows the modification induced by stellar toroidal oscillations on the accelerating 
electric field and on the charge density. We also find that, after decomposing the oscillation 
velocity in terms of spherical harmonics, the first few modes with m = 0, 1 are responsible 
for energy losses that are almost linearly dependent on the amplitude of the oscillation and 
that, for the mode (/, m) = (2, 1), can be a factor ~ 8 larger than the rotational energy losses, 
even for a velocity oscillation amplitude at the star surface as small as r\ = 0.05 O R. The 
results obtained in this paper clarify the extent to which stellar oscillations are reflected in the 
time variation of the physical properties at the surface of the rotating neutron star, mainly by 
showing the existence of a relation between PP and the oscillation amplitude. Finally, we 
propose a qualitative model for the explanation of the phenomenology of intermittent pulsars 
in terms of stellar oscillations that are periodically excited by star glitches. 

Key words: MHD: pulsars — general — relativity — oscillations — stars: neutron — plasma 
magnetosphere 



1 INTRODUCTION 

The theoretical study of radio pulsars dates back to the work of Goldreich & Julian 1 19691) who first suggested the existence of a magneto- 
sphere with a charge-separated plasma around rotating magnetized neutron stars. A spinning magnetized neutron star generates huge potential 
di fferences betwee n di fferent parts of its surface. The ca scade generation of electron-positron plasma in the polar cap region already proposed 
by Sturrockl 1 1971 ) and Ruderman & Sutherland 1 1975 ) requires that the magnetosphere of a neutron star is filled with plasma, thus screening 
the longitudinal electric field and bringing the plasma into co-rotation with the neutron star. Because co-rotation is not possible outside the 
light cylinder (the radius Rlc = c/Q at which the co-rotation speed equals the speed of light), essentially two different topologies of the 
magnetic field lines are naturally produced: closed lines, namely those returning to the stellar surface, and open lines, i.e. those crossing the 
light cylinder and going to infinity. As a result, plasma may leave the neutron star along the open field lines and it is generally thought that 
pulsar radio emission is produced in the region of open field lines well inside the light cylinder and within a g i ven ang l e from the polar axis. 



ar radio emission is produced in the region ot open field lines well inside tne lignt cylinder and witnin a given ang] 
Beside the seminal papers by Goldreich & Julian! 1 1969h . Sturrock ( 1971 ), Ruderman & Sutherland 1 1975 ), Mestell 



1 197 if) and lArons & S charlemann 



( 1979J), pulsa r magnetospher es h ave been investigated by a la rge number of authors over the years. We only mention here the reviews by 
lArons! Jl99lh . Mestell 1 1992 ) and Muslimov & Hardingl 1 1997 ). where subsequent achievements and so me new ideas h ave been presented. 
Thorough description of known magnetosphere properties may be found, for example, in the book of iBeskin J2009h . It should also be 
mentioned that in the last few years time dependent numerical simulations of neutron star magnetospheres have been proposed as a new 
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promising tool for investigating the complex physics of these systems. At least qualitatively, the numeric al approach has confirmed the most 
funda mental features of what was expected from the stationary solution of the Gra d-Shavranov equation I Contopoulos et al.lll999l ; Gruzinov 
2005 ), such as the existence of closed magnetic field lines up to the light cylinder l Komissarovlbood . McKinnevlbood) . or the scaling of the 



spin dow n luminosity on th e angular velocity and on the inclination angle of the neutron star angular momentum with respect to its magnetic 
moment I Spifkovskv 2006h . In spite of this spectacular progress, however, the numerical approach still suffers from some serious limitations, 
such as the lack of a unified scheme in which both the force free regime and the plasma regime of magnetohydrodynamics are simultaneously 
taken into account, or the lack of a consistent treatment of resistive effects in the current sheet. 

The analytic approach, on the other hand, can still provide a deep understanding of pulsar physics. In parti cular, a lot of a ttention 
has been paid to the ex istence of a strong electric field induced by the rotation of the star, as already noticed by Deutschl 1 1955). More 
recently, Beskinl 1 1990h and, independently, Muslimov & Tsygan 1 1990l) were the first to find that the frame dragging induced by general 
relativistic effects provides a source of additional electric field contributing to particle acceleration in the polar cap region. The accelerating 
component (parallel to the magnetic field) of the electric field is driven by deviations of the space densit y charge from the Goldreich- Julian 
(GJ) charge density, which is determined by the magnetic field geometry. As noted by several autho rs i Beskin 1990l : Muslimov & Tsygan 
199dlMuslimov & Harding||l997l : lDvks et al]|200ll ; lMofiz & Ahm edov 20o"(il: lMorozova et alj|2008h . the corrections of general relativity in 



the plasma magnetosphere of rotating neutron stars are first-order in the angular velocity of the dragging of inertial frames and have to be 
carefully included in any self-consistent model of pulsar magnetosphere, especially when computing the resulting electromagnetic radiation. 

Tightly related to this aspect is the possibility that neutron star oscillations, most likely excited during a glitch phenomenon (sudden 
change of the rotational period), propagate into the magnetosphere, thus affecting the acceleration properties in the polar cap region . The 
first attempt to generalize the Goldreich- Julian formalism to the case of an oscillating neutron star was made bv lTimokhin et al. I hood) , who 
developed a general procedure for calculating the GJ charge density in the near zone of an oscillating neutron star. Using this procedure, the 
GJ charge density and the electromagnetic energy losses were computed for the case of toroidal osc i llation s at the neutro n star surface. A 
simila r approach has been recently extended to the general relativistic context by Abdikamalov et al. 1 20091) who, just like Timokhin et alj 
i 200Ch . based their results on the so called low current density approximation, i.e. on the assumption that the magnetic field is mainly 
produced by volume currents inside the neutron star and by surface currents on its surface, while the magnetic field due to magnetospheric 
currents can be neglected. In the paper of lAbdikamalov et al.1 Jiooj) the influence of oscillations to the magnetos phere electrodynamics was 
considered for the case of a non-rotating Schwarzschild star. In the present paper we apply some of the results of lAbdikamalov etaL 1 20091) 
to investigate how oscillations, produced at the star surface, refle ct in the energy losses from the polar cap region of the magnetosphere of 
slowly rotating neutron star. In this respect we extend the work of lMuslimov & Hardingl d 19971) by performing a local analysis in the domain 
of open magnetic field lines in the inner magnetosphere and taking into account the effects of toroidal oscillations excited at the star surface. 

The plan of the paper is as follows. In Sec. [2] we provide the minimum general relativistic formalism for understanding neutron star 
electrodynamics and we perform a detailed analysis of the GJ charge density of slowly rotating and oscillating neutron star. In Sec.[5]we derive 
a version of the Poisson equation that takes into account both general relativistic effects and the oscillating behavior of the magnetosphere of 
the rotating star. Sec. [4] on the other hand, is devoted to the computation of the energy losses induced by oscillations together with rotation. 
In Sec.[5]we propose and motivate a suggestive idea to explain the phenomenology of intermittent pulsars in terms of the excitation of stellar 
oscillations. Finally, Sec.|6]contains the conclusions of our work. 

Throughout, we assume a signature { — ,+,+,+} for the space-time metric and we use Greek letters (running from to 3) for four- 
dimensional space-time tensor components, while Latin letters (running from 1 to 3) will be employed for three-dimensional spatial tensor 
components. Moreover, we set c = G = 1 (however, for those expressions with an astrophysical application we have written the speed of 
light explicitly). 



2 GOLDREICH-JULIAN RELATIVISTIC CHARGE DENSITY 



In the slow limit approximation, the spacetime around a rota ting neutron star of total mass M, angular momentum J and angular velocity O 
is given by faartle & Thornelll968l : Eandau & Lifshitzil2004l) 



-N 2 dt 2 + N~ 2 dr 2 + r 2 



9 2 + sin 2 



2oji.Tr 2 sin 2 9d<f)dt 



(1) 



where N = (l-2Af/r) 1/2 is the gravitational lapse function, while ljlt = 2aM/r 3 is the Lense-Thirring angular velocity, which represents 
the angular velocity of a freely falling inertial frame. The specific angular momentum a, on the other hand, is defined as a — J/M, J is 
the an gular momentum of the star. We note that the metric is split according to the 3+1 formalism of general relativity i Arnowitt et alj 
1962h . which admits a natural Eulerian observer, also called the ZAMO (zero angular momentum observer), with four-velocity n a given by 



n a = {-N, 0,0,0} 



(2) 



In the rest of our discussion, when we introduce any three- vector, like for instance the electric field, we assume that it is defined in the locally 
flat spacetime of the ZAMO observer, and we denote its orthonormal components with hat superscripts. 

From the system of Maxwell equations and after assuming the magnetic field of a neutron star to be stationary in the co-rotating frame, 
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Muslimov & Tsvgar] 1 1992 ) derived the following Poisson equation for the scalar potential $ 



V-|-V$ 



-47r(p - pgj) 



(3) 



where p — pc,i is the effective space charge density responsible for the generation of an unscreened el ectric field parallel to the m agnetic 
field, while pgj is the Goldreich- Julian charge density that we discuss below. In their pioneering work. lGoldreich & Julian! J1969T) showed 
that a strongly magnetized and highly conducting neutron star, rotating about the magnetic axis, would spontaneously build up a charged 
magnetosphere. In a nutshell, the argument is the following: if a magnetized rotating neutron star is placed in vacuum, enormous unbalanced 
electric forces parallel to the magnetic field B would set up at the surface of the star, e xtract ing charges from the surface into the external 
vacuum region, thus producing a filled magnetosphere. Therefore, Goldreich & Julian 1 19691 ) hypothesized that a far better approximation 
for the magnetosphere would be obtained by shorting out the component of the electric field E along B. The m agnetospheric charges tha t 
maintain E ■ B = are themselves subject to the E x B drift that sets them into co-rotation with the star jMofiz & AhmedovfeoOol) . 
A derivation of the Goldreich- Julian charge density in the presence of oscillations but in the Newtonian framework has been performed 
by Timokhir] 1 2007 ). Here we discuss the corrections to the standard Goldreich- Julian charge density when both relativistic effects and stellar 
oscillations are taken into account. Our starting point is the general expression for the Goldreich- Julian charge density that takes into account 
the contribution of the electric field induced by arbitrary stellar oscillations, i.e., 



Pgj 



47TC 



— (u — w) x B + j^Sv x B 



47TC 



N V f 3 / 



u x B + j^Sv x B 



(4) 



where u — w = (fl — cjlt) r sin 9e$, while Sv is the oscillation velocity. Moreover, R is the star radius, r = r/R is the dimensionless radial 
coordinate, e = 2M/R is the compactness parameter, p = I/Iq is the moment of inertia of the star in units of I — MR 2 and n = e/3. 



W e apply expression l(4} to the case of toroidal oscillations, whose velocity has components [see, for example, Eq. (13.71) of lUnno et al 



8v l = 



0, ^—^d,pYi, m ,( 
sm 9 



),-d 6 Y llm ,(B,<l>)\f,{r)e 



(5) 



where u is the real part of the oscillation frequency, while i) is the radial eigenfunction expressing the amplitude of the oscillation. We have 
used multipolar indices I' and m' to distinguish the harmonic dependence of the velocity perturbations from the harmonic dependence, in 
terms of I and m, of the electromagnetic fields, since these indices are in general distinct. As usual, the spherical orthonormal functions 
Yim(0, <j>) are the eigenfunctions of the Laplacian in spherical coordinates. They are given by 

1 



Yu 



where the functions 6>; m (cos 9) satisfy the differential equation 



d 



dQ [m (cos9) 
d9 



m 
sin 2 6> 



°e im (cos(9) , 



e im (cos 9) + 1(1 + l)e !m (cos 9) = 



and can be written as 



e im (cos9) = (-i) r 



21 + 1(1- m)\ 



Pr (cos 9), 



(6) 



(7) 



(8) 



(I + m) 

where P ; m (cos 9) are the Legendre polynomials. For simplicity, we limit our attention to the case in which the magnetic moment /i of the star 
is aligned with its angular momentum, and, furthermore, we assum e that the magnetic field is a dipolar one, with orthonormal components 
that take the form ( Muslimov & Tsyganlll99ilGinzburg V. £ll964l) 

3 



nf D J\ T ) —3 o r,0 

B =B j^-r cos 9 , B 



/(I) + (1 - e/f)f(l) 



' sin 9 . 



(9) 



(10) 



where 

and Bo = 2/i/i? 3 is the Newtonian value of the magnetic field at the pole of the star. When the oscillation velocity and the magnetic field of 
the star are given, respectively, by Eq. 10 and Eq. l|9}, then Eq. (|4j provides the following modified Goldreich- Julian charge density 

SIB 1 f(r) 



PGJ 



2ttc Nf s f(l) 

1 1 R 

^TcR^ B ° e 



1-4 



Nf(l) 



9 7^ sin ^ yi '-'( e '^ 



;Y Vm ,(9,4>) 



+ 



N 
~2 



-2 



f(f) 3 
/(l) + (1 - e/f)f(l) 



df N 



(11) 



From the point of view of the emission of energy, that we consider with greater detail in Sec. [4] the most interesting region of the magneto- 
sphere is the so called polar cap region, i.e. the region where the magnetic field lines remain open and at distances from the star surface much 
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smaller than the light cylinder radius. If we denote by Oo (computed as explained in Sec. [4} the co-latitude of the last closed magnetic field 
line a t the star surface, then the angle 6 of the last closed magnetic field line as a function of f is well approximated as (M uslimov & Tsygan 
1992D 

1/2 



6(f) 



rim 



sin 6o 



(12) 



It should be remarked that on the surface of the star, and even within reasonably large distances far from it, the aperture angles of open 
magnetic field lines remain small. For example, for a neutron star with e = 0.3 and Go = 0.087 « 5°, we find G = 0.197 ~ 11.3° at 
r = 4R. As a result, in this approximation the last term in the curl brackets on the right hand side of Eq. J 1 1 b can be neglected, since it 
contains 9, with 9 ^ G(f), to a second power larger than other terms in the same brackets. Taking into account Eq. ((6} and 10, in the limit 
of small angles 6, we obtain from i ll It the following equation for the Goldreich-Iulian charge density 



PGJ = PG.L0 + 5pG3, I'm' = 



VB 1 /(f) 



1 - 



1 1 B e~ 



i m 



V(f)l'(l' + 1)Y V 



(13) 



2ttc Nf 3 /(l) V f 3 ) 4-KcRf 4 G 2 (f) N f(l) 
where pgj,o is the Goldreich- Julian charge density of a slowly rotating neutron star while <5pgj, I'm' is the correction induced by oscillations 
We are here interested in analyzing the Goldreich-Julian charge density of the first few modes, namely those with (?', m') given by (0, 0) 
(1, 0), (1, 1), (2, 0) and (2, 1). To this extent, however, we greatly simplify our calculations by approximating Vj/m' ifii 4 1 ) ~ ^-Vm 
where the terms Ai m (<f)) have real parts given by 



An 



1 



A m = \ -r- 



/47T ' V 47T 

From d!3t we can compute the ratio 



5p 



GJ 1' 



'/PGJ,0 = — 



K 



2f2 — m/2 



U 



m 



A 2 i = -3 



+ l)A llm ,{4>) 



24tt 



■jk) 



(14) 



(15) 



where we have posed 77(f) m fj(l), which amounts to the assumption that the oscillation amplitude maintains the value it has on the surface of 
the star, at least within small distances far from it, as we are considering here. Moreover, we have introduced the small number K — fj(l)/QR 
to parametrize the amplitude of the oscillation. Finally, have considered G ~ 6. One can easily see that <5pcj I'm' = f° r me mode (0, 0). 
FigureQ] on the other hand, shows the ratios SpQj y m i / pG3,a for the other four modes, computed at t = 0. When plotting these graphs we 
have used the following typical set of parameters: k = 0.15, e = 1/3, K = 0.01, Go = 0.008, Q. = lrad/iQ As it is clear from Fig.Q] the 
oscillation induced Goldreich-Julian charge density not only can be a significant part of pgj.o, but can even prevail several hundred times 
over it, for example, for the mode (2, 0). An exception is represented by the mode (1, 1), for which \5pgj y m i \ < |pgj,o| even very close to 
the star. The influence of oscillations is greater near the surface of the star, which is the most interesting region of the magnetosphere, while 
it decreases far from it. On the other hand, the relativistic charge density p that enters Eq. Q is proportional to the inten sity of the magnetic 
field through a proportionality coefficient that is constant along the given magnetic field line JlVIuslimov & Tsyganlll99ll) , i.e. 

flB 1 /(f) 



P = 



+ e 



(16) 



2ttc Nf 3 /(l) 

where we have introduced the variable £ = #/G(f), and where A(£) and a(f, <j>) are unknown functions to be speci fied from the boundary 
condit ions. The computation of A(^), which corresponds to the case of pure rotation, has already been performed bv lMuslimov & Harding 
1 19970 (see their Eq. (58)), showing that A(£) w k — 1. The computation of a(£, (f>), on the other hand, is discussed in Sec. l3.1l below. 



3 POISSON EQUATION 

After inserting d 1 3 1 > and \\6\ into Eq. y), we find the following expression for the Poisson equation 

1 



r 2 dr V dr 



+ 



Nf 2 9 



89 J 6 dcf> 2 



-47T 



SIBo 1 /(f) 



3 V® 



+ e~ 



f a(e,0) 



4ty 



$ = 

fiBo 1 /(f) 



1) V f 3 J 



1 1 B e~ 



1 m f,{f)l'{l' + l)Y v 



2nc Nf 3 /(l) L vs; vs ' ry J 2ttc Nf 3 /(l) V f 3 1 c Rf A G 2 (f) N /(l) 

Solving such equation represents a complicated task, and our strategy consists in, first of all, expanding the scalar potential $ as 

fc(t, f , £, 4>) = Mr, + e-^S^f, £, 4>) , 



(17) 



(18) 



We will show in Sec[4]below that K and ©o are not independent of each other and that 0o does also depend on the indices (' and m' . However, we have 
chosen to use a single value 0o = 0.008 for all of the plots reported in Fig.[TJas this nevertheless represents a mean value typical of a standard astrophysical 
situation. 
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Figure 1. Ratio SpQj i' m '/pGJ,0 f° r the mode (1, 0) (left top panel), for the mode (1, 1) (left bottom panel), for the mode (2, 0) (right top panel) and for 
the mode (2, 1) (right bottom panel). The representative star parameters are chosen to be k = 0.15, e = 1/3, K = 0.01, 0q = 0.008, Q = lrad/s. 



where the first term 5>o(f , £) corresponds to the case of a non-oscillating but rota ting star. In this case, in fact, we know that $o must satisfy 
the following equation [see, for example, Eq. (37) of lMuslimov & Tsvga 



R~ 



1 d 



r 2 dr 



d 



dr 







Nf 2 6 86 V 06 



'4 1 



2QB /(f) 



/(I) \N 



^(0 + ^(1-- 



(19) 



The next step consists in expanding in terms of spherical harmonics both the second term of J 1 3b . representing the perturbation to the 
non-oscillating case, and the function 5(£, cj)), which contains the <j> dependence of the charge density d!6t . Namely, we write 



1=0 m=-l 
oc I 



(20) 
(21) 



At this point, lengthy but straightforward calculations allow to derive from Eq. d!7t the following equation for the perturbation 5$(r, ( 



S$, m Y, m = -4tt 



ftBp 1 /(f) 
2-kc JVf 3 /(l) 



1 1 



Bn 



1 /(f) -j7(f)/'(r + i)y ! ,. 



(22) 



ci?f 4 e 2 (f) jv/(i) 

We now introduce 5Fi m (r) = f5$; m (f), and we exploit the fact that the spherical harmonics Yi m form a set of orthogonal basis functions. 
This means that in Eq. \22\ we must first fix /' = I and m! = m and we then equal the coefficients of each basis function Yi m on the left and 
on the right hand side. In this way we finally obtain the Poisson equation in the form 



R- 



2 N_ ( jf 



1(1 + 1) 



\df 2 N 2 f 2 e 2 (f) 



SFi m {f) 



Bo /(f) 



1 fj(r)l(l + l) 
Rf 9 2 (f) 



(23) 



cNf s /(l) 

It is worth stressing that Eq. \23\ is valid for small polar angles 9, but for any distance f from the star surface. The limit of Eq. d23 b for small 
f is considered in the rest of this Section. 
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3.1 Solution close to the star surface 



As a first example, we wish to compute the solution of the Poisson equation close to the star surface, where z = r- l<l, while imposing 
no restrictions on the aperture angle Bo of the last closed magnetic field line. In this case Eq. i23\ for the unknown SFi m becomes 



d 2 _ Z(Z + 1) 
dz 2 ~ (1 - e)e ( 2 , 



x „ , . 2Q,B Q R 2 1-2z. 
d±>im{z) = aw, 



BpR 1 - 3z 
cQ 2 1-e 



mm +i) 



(24) 



which represents the extension of Eq. (44) of Muslimov & Tsyganl 1 1992 ) to account for the presence of oscillations. The general solution of 
Eq. d24b . which is an inhomogeneous differential equation, is 



5F lm = Ce 



VHi + i) 



B R 



1(1 + 1) c 



2flR(l - 2z)a lm + l JL^±rjQ.)(l - 3z) 



(25) 



The constant C and the coefficients a; m may be found after imposing physically motivated boundary conditions. In particular, at the star 
surface we require both equipotentiality and absence of a steady state electric field, which amounts to the two conditions 

dSFlm 



SF,, 



z =o = . 



dz 



From them, simple calculations allow to derive 



dim = 



!(/ + !) 



1 - 



= 0. 



c = 



B a R 



^/l(l + l)-2Q VT=l_ 

The final solution for the scalar potential near the surface of the oscillating rotating neutron star has therefore the form 

oo I 



(26) 



(27) 



*oC.fl+e-*-J^(l)£ E 



QoVT 



1 = m=-l 



^l(l + l)-2Q y/T=< 



■e vr=W ' + 1- 



and the corresponding expression for the accelerating component of the electric field is 

oo I 



En = Ei 



^ e » ( J -lj>Yi m (£,< 



(28) 



(29) 



where 



is the accelerating field that is present even in the absence of oscillations. 



(30) 



3.2 Solution in the polar cap region 



The solution in the polar cap region, namely when 6o < r - 1< Rlc/R, can be obtained from Eq. j28) and Eq. l |29t in the limit of small 
6 , or, directly, from Eq. (23]l. In this case, in fact, \d 2 SF im /df 2 \ < 1(1 + l)\SFi m \/N 2 f 2 B 2 (f), and Eq. {23) therefore reduces to 



Z(7 + l) 



2QB 



fir)- ,1 So 
-aim + 



fir) 



*j(f)i(l + l) 



R 2 Nf 3 Q 2 (f) c Nf 3 /(l) c9 2 (f) i?iVf 4 /(l) 

After using O as given by Eq. dl2t and with the same coefficients 2; m expressed by Eq. l !27t . we can obtain 



(31) 



SFlm = 



B R 



-rfftl) — ;====: . h 



y/l(l + 1) - 2vT^7e„ /(I) r 
This allows to derive both the electric potential and the accelerating component of the electric field as! 



*o(r,£) 



-Bo-R 



fir) vir) 



/(I) 



E E 



(32) 



(33) 



Eo-^— ^ /(f) ^ 



E E ( 34 ) 

i=0 m=— ! 

It is interesting to compare the seco nd term on the right hand si de of Eq. ( 134b . namely the contribution due to the oscillations, with Eo as 
reported, for example, in Eq. (55) of iMuslimov & Tsyganl ( 119921) . We find 



5E 



E 



lm —itot 

= e 



2 *7(1) 
3 Q.Rn 



d_ //(f) 1 

df v/(i) f 2 



rOo 



■^imO) 



(35) 



For small values of 0o we use the approximation + 1) - 3 % /T^70 o )/( + 1) - 2 % /T^7eo) ~ 1. 
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Figure 2. Ratio 5E i m /Eo for the mode (1, 0) (left top panel), for the mode (1, 1) (left bottom panel), for the mode (2, 0) (right top panel) and for the mode 
(2, 1) (right bottom panel). 



that we plot in Fig. [2] for the different modes considered so far and computed at t = using the same set of parameters as for the plots in 
FigQ] Moreover, we have posed £ = 1 /2, i.e. we have considered the middle magnetic field line between the polar axis and the edge of the 
polar cap. The plots demonstrate the significant influence of oscillations to the electric field of a pulsar. In particular, because Bo is supposed 
to be small, the ratio 8E i m /Eo can be very large for modes with m < 2. Indeed, the modes (I, m) — (0, 0) , (2, 0) induce an electric field in 
the opposite direction to Eo, which is three orders of magnitude larger in modulus than Eo, In general, the ratio \ 8E i m /Eo\ increases when 
increasing the distance from the star surface. 

In a recent work aimed at quantifying the impact of neutron star oscillations on the accelerating electric field, Timokhir] ( 2007 ) consid- 
ered the whole range of angles 6 and more complicated cases of spheroidal oscillations. He found that the contribution of oscillations to such 
electric field may be either positive or negative and that this contribution is substantial for values of / of several hundreds. Indeed, as shown 
in the Sec.|4]below, the appearance of the small parameter Oo in the computation of the energy losses makes astrophysically relevant even 
the modes with small values of I and m. 



4 ENERGY LOSSES 



In this Sect ion we calculate the energ y losses from the polar cup region of a rotating and oscillating neutron star, using several results 

(2009) for the spherical Schwarzschild star. The tota l energy loss from the open field lines region, averaged 



Abdikamalov et al 



presented m[ 

over an oscillation period and carried out by the outflowing plasma, is determined as ( Timokhin et al.ll2~000h 



dt 



d8\jf m (R, 6, <t>)Ae lrn (8, <f>)\R 2 sin( 



where Ae(9, <j>) is the work done by the electric field to move a unit charge to the point with coordinate (R, t 

f0 



Aei m = RNi 



Eg j im(R,6' ',<!>)<& , 



(36) 



(37) 



while jf m is the electric current density and r is the period of oscillations. It should be noted that despite some terms in the expression for 
the energy losses will contain the factor e~ , which mathematically gives a zero net result when averaged over time, physically these terms 
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will give a non zero net contribution. That is because the average crossing time of the acceleration zone by particles which are responsible 
for the energy losses is much shorter than the oscillation period (typical frequencies of oscillations are of the order of 10 — 50kHz). In 
practice, therefore, the particles do not return back to the star (see also the discussion of Timokhin et al, 1 2000l)). If the po l ar cap region is 



in the conditions of complete charge separation, the current density is well approximated as i Ruderman & Sutherland 1975 ; Timokhin et al 
2000h 



f =± Pgjc . (38) 

In order to compute the integral in B6t . we need to provide an estimate of the angle Oo of the last closed magnetic field line at the surface 
of the star. After using R a to denote the radial coordinate of the point where the last closed magnetic field line crosses the equatorial plane 
(note that R a 3> R), we derive Bo from the condition 

e p i (R a , tt/2, <j>) = e em (R a , it/2, <f>) , (39) 

where e p i and e em are the kinetic energy density of the outflowing plasma and the energy density of the magnetic field, respectively. The un- 
derlining idea, in fact, is that at the last closed magnetic field line equipartition of energy exists between the magnetic field and the plasma, and 



for co nvenience, the condition is evaluated at the equator. The two energies e p i and e em have already been computed bv lAbdikamalov et al 
(2009) [see their Eq. (100) and (101)] and are 



£pl(Ra, 7T/2, 



i(Ra,n/2, 



_J_Nr_R_ 

2/(1) N Ra R a 

N Ra R 6 2 
32tt R n ' 



Aef 



(40) 



(41) 



where N R = \Jl- 2M/R, N Ra = \J\ - 2M/R a w 1. The angle 9 does not explicitly appear in (gojl and (HQ, but rather implicitly. In 
fact, for a dipole magnetic field f(r) sin 2 6/r = const and therefore 

R /(i) 

Ra 



(42) 



f(Ra/R) 

where f(R a /R) is close to unity. In order to proceed with the computation of Oo, and hence of the energy loss, we need the electric field 
component 



(E°) 



GJ Ir, 



Nc 
1 



u + 8v\ B, 



— [siRf sine (l- 21) - d s Yi m (e,<t>)fi(r)\ B 



/(i) ^ 3 



from which we can compute Ae ; m through d37t . This allows us to obtain (for small angles 9) 



Aei 



RN R B 



rn r b e 2 



QR(1 - 
[flR(l 



A lm 6 m rj(l) 



k) - 2A J0 ij(l)] 



for m/0 



for m = . 



(43) 



(44) 



(45) 



c 2 

Taking into account equations l |39l l, d40t and J4U . we can obtain the expression for the angle Oo, which is an algebraic equation in the case 
m/0 while it has a closed form in the case m — 0, i.e. 



RN R 



QR(1 - k) 



A lm O™rj(l] 



n(i 



c 2cR 



8 



/ 4 (i)e§ 



©0 



2 VR^NK 



[Q.R(l -k)- 2A lo fj(l)] 



1/4 



^ c (l-K) + j- c ^fj(l)l(l + l)A l0 



1/4 



for m/0, 
for m = 



When ?n = and O = 0, hence in the axisymmetric case with no rotation, we get the expression 



6o I 



4A* i(Z + l)JW(l) 



/ 4 (1) 



This estimate is slightly different from that obtained by Abdikamalov et al. 1 20091) 



6o = 2N R < 



1/4 



V(1)A% 



/(I) 



1/2 



(46) 
(47) 

(48) 

(49) 



firstly because of the different expansion of the spherical harmonics, which in Abdikamalov et al. i2009j) was chosen to be Yi m (9, 4>) w 
A \]l (4>)Q m + A^ (4>)8 m+2 , and, secondly, because of the different determination of Bo (equations (73)-(75) for magnetic field in the paper 
of Abdikamalov et al. | j2009h do not contain /(l) in the denominator). On the other hand, when ra — and rj — 0, hence in the axisymmetric 
case with no oscillations, we get the expression, 



0o|m=O,77=O 



4N R Q 2 R 2 
/ 4 (1) ^ ' 



(50) 
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which represents the correct general relativistic extension of the expression reported bv lMuslimov & Tsvganl Jl992h . Inserting the angles Oo 
defined by the equations d46b and d47b into d36b . we derive the total energy losses from the polar cap region for a rotating and oscillating 
magnetized neutron star as 



r 3 n r bZ 



2tt 



SI 



1 O m +4 



1 1 

2cRN R 



A 2 m (0)tf (1)1(1 + 1) 



(-) 



2m + 2 



2m + 2 



, R3 ^ B2 ° ^ m{1 2^0(^)77(1)] -k) + Yc^ mi + 1)AloW 



for m # , (51) 



for m = . (52) 



Equations d51b and d52b just derived are astrophysically very relevant and deserve some comments. In the first place, it is interesting to note 
that d52b takes a simpler form if only the linear terms in the amplitude of the stellar oscillation are retained. In this case, in fact, we find 

fj(l) Aw I 2 + 1-2] 



R 4 B 2 Sl 2 
8c ' 



k) e 



i + 



SIR 1 



(53) 



We can highlight the corrections with respect to the energy losses in the Newtonian case and in the absence of oscillations if we replace Oo 
in d53b with one of the expressions computed above. For simplicity, we consider the case given by d50t . namely the case of pure rotation, and 
we obtain 



3(1 



N R 
/ 4 (1) 



1 + 



77(1) A l0 I 2 + 1-2 



SIR 1 — re 



(E ro 



tJNewt j 



where (E ro t 



is the standard Newtonian expression for the magneto-dipole losses in flat space-time approximation 

■ _ 1 Sl 4 B 2 R 6 

{ Cjrot) Newt — — o ■ 

6 cr 



(54) 



(55) 



The astrophysical relevance of Eq. d54b becomes even more transparent when it is rewritten in terms of the pulsar observables P, i.e. the 
period, and P = dP/dt. To this extent we first recall two standard relations of pulsar physics. The first one is the relation between luminosity 
and spin-down, namely 

L = -ISISI (56) 

where / is the general relativistic moment of inertia of the star ( Rezzolla & J. Ah medov 2004 ). The second (Newtonian) relation is between 
the spin down and the intensity of the magnetic field, i.e. 

(2t^ \ R 6 B 2 



(PPh 



V 3c 3 



where / is the classical moment of inertia of the star. From these two relations and from Eq. ( 154b we deduce 



(PP) 



1 + 



jj(l) Aio I 2 + 1-2 



SIR 1 



/ N E 



r/ 4 (i) 



{pp)i 



(57) 



(58) 



Pulsar periods P and spin-down rates P are very precisely measured quantiti es for a large number of pulsars. iKaspi et al.l d2004f) , for instance, 
report &P — P diagram of 1403 cataloged rotation-powered pulsars (see also lAronsI d2007h ). Thus, expression d58t for PP can in principle be 
applied to existing observations for detecting the possible effects of stellar surface oscillations. The main difficulty encountered in this kind 
of analysis rests in the low accuracy measurements of the stellar radius, which reflects in poor estimates of the moment of inertia. However, 
the predictive potentiality of d58b is unquestionable, and once the moment of inertia of neuron stars is determined more precisely, it will 
provide detailed information about stellar oscillations starting from fundamental observational properties about pulsar timing. 
We now discuss the energy losses in the case m 7^ 0, and to this extent we rewrite Eq. d51b as 



— L r ot 



dS— 

Y 2n 



l + K 



1 - K 



l(l + l)A lm {<t>) 



-K- 



8 



1 - K 



A, 1^ 



R 4 B 2 Sl 2 
8c 



- k) 2 Q 4 I d0— 



K 2 ti(l + l)Af m (0) 



<-);, 



2m 



where Lrot denotes pure rotational energy losses 



R 4 B 2 Sl 2 . ,2 Q 4 
— 1 1 — K J "0 



(59) 



(60) 



We are here interested again in analyzing the energy losses of the first few modes, namely those with (I, m) given by (0, 0), (1, 1), (2, 0) and 
(2, 1). Eq. d59b contains two major contributions. The first contribution includes the rotational energy losses plus terms linear in K which 
are not present neither in pure rotational nor in pure oscillatory regime. The second contribution, on the other hand, is the only one present in 
the pure oscillatory regime. More specifically, in the pure rotation regime (7) = 0, SI 7^ 0), only the first term in the curl brackets of Eq. d59b 
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1-2, m-1 
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Figure 3. The angle (in deg units) of the last open magnetic field line ©o as a function of the parameter K = ij(l)/QR for modes (1,1) (continuous red 
line) and (2, 1) (dotted blue line). The representative parameters of the star are R = 10km, Q = lrad/s, e = 1/3. 



survives, and we simply recover L — L rot . In the opposite regime of pure oscillation (77 ^ 0, £1 = 0) the whole first contribution to the right 
hand side of Eq. i59l vanishes and only the second contribution proportional to K 2 survives. Finally, in the mixed regime with both fj 7^ 
and £1 7- 0, all of the terms must be included. 

In order to appreciate the dependence of the energy losses on the oscillation amplitudes, we have computed the ratio L m ^o /L rot as a 
function of K for all of the modes mentioned above. As typical and representative parameters of the rotating neutron star we have chosen 
R = 10km, Q = lrad/s, e — 1/3. We first computed the angle Bo after applying a standard root solver to Eq. d46t . Graphs of Qq as a 
function of the parameter K are presented in Fig.[5]for the modes (1, 1) and (2, 1) and they show that the size of polar cap increases with the 
amplitude of stellar oscillations. Fig. [4] on the other hand, reports the ratio L m ^o/ L ro t. Interestingly, it follows from Eq. d59b that modes 
with m = 1 have the small parameter Oo to a negative power. Therefore, for such modes the energy losses due to oscillations may exceed 
significantly the energy losses due to pure rotation, even for relatively small K. This effect is indeed shown in the left panel of Fig. [4] which 
reports the ratio L m -t / L ro ± for the modes (1,1) and (2, 1). The right panel of Fig. [4] on the other hand, reports the ratio L m ^IL rot for the 
modes (0, 0) and (2, (ijl and they turn out to be a factor ten smaller than the ones reported in the left panel, for the reason explained above. 
We note, on the contrary, that the energy losses due to oscillations are smaller than those due to rotation for the mode (0, 0). Modes with 
m > 1 do not contain Oo to a negative power. Although modes with higher m practically do not give any contribution to the energy losses, 
we should remark that for modes with m > 3 the angle Oo may not be small, thus requiring an alternative approach to the one presented in 
this paper. Since the energy losses of the pulsar are proportional to the spin down rate, as it is clear from Eq. d56l l. the same graphs obtained 
for the ratios L m / L rot describe also the ratio Cli m /Cl rot , where fi; m denotes the time derivative of the pulsar rotation frequency when the 
star oscillates with the mode (I, m) while tl ro t corresponds to the case of pure rotation. Thus, it follows from Fig.|4]that an observable value 
of Qi m /Q, ro t ~ 1.5 may be reached for a value of K as small as K ~ 0.03 for the mode (1, 1) and for K ~ 0.01 for the mode (2, 1). 

Finally, it is interesting to note that the comparison of our results with those of iTimokhin et al. I hood) can be done only with caution, 
since he considered the whole range of angles 9, more complicated cases of spheroidal oscillations, and he did not show the ratio L m -t / L ro ±, 
but rather the ratio bet ween plasma energy losses and vacuum energy losses, both in the presence of oscillation^- It should also be noted that 
20001) performed estimations rather than exact calculations. 



Timokhin et al 



5 CONNECTION TO THE PHENOMENOLOGY OF PART-TIME PULSARS 

Recently, Rea et al] d2008h and lLvnel d2009h reported that a previously known pulsar, PSR B1931+24, with a spin period of 813 ms at the 
relatively large distance of ~ 4.6kpc, when monitored long enough, shows an intermittent radio emission, consisting of an active ON state, 
lasting for 5 — 10 d, followed by a sharp transition (happening in less than 10 s) to an OFF state during which the pulsar remains undetectable 
for 25 — 35 d. More interestingly, the spin-down rate during the ON state is !>on = —16.3(4) x 10 _15 Hz s _1 , while, when measured over 
longer periods, the average value of the spin down is sensibly smaller. This is compatible with a picture in which the spin down rate during 



3 The mode (1, 0) does not have linear contributions in K to the energy losses, as it is follows from (54). 

4 What he found is that for some modes, like for instance (1, 0) or (I ^ 2, m = 0, 1) the plasma energy losses are larger than the vacuum ones even for 
small oscillation amplitudes. 
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Figure 4. Left Panel: The ratio L m ^o/ L ro t as a function of parameter K = r){l) /QR for modes (1,1) (continuous red line) and (2, 1) (dotted blue line). 
Right Panel: The ratio L m —o/L ro t as a function of parameter K = fj(l)/QR for modes (0, 0) (continuous red line) and (2, 0) (dotted blue line). 



the OFF state is different from that in the ON state, and, in particular, it is a 50% factor smaller, namely !>off = —10.8(2) x 10 15 Hz s . 
A further search for similar objects from the Park es Multi-Beam Survey data revealed at least four additional objects presenting properties 
similar to those of PSR B1931+24 teeckejEooi) . like, for instance, PSR J1832+0031 with an ON state of ~ 300 d and an OFF state of 
~ 700 d. Understanding the physical mechanism responsible for such a remarkable phenomenology as well as the relationship between 
intermittent (or "part-time") pulsars and conventional radio pulsars is of course of great interest, especially as it may help clarifying aspects 
that are still obscure about pulsar radio emission. 

All of the (few) models presented so far for explaining the phenomenology of intermittent radio pulsars are based on the common idea 
that intermittent pulsars are isolated neutron stars, simi lar to conventi onal radio pulsars. One of the first ideas was that this effect could be 
similar to nulling, already reported several years ago bv lBackenJl97(il) . However, the nulling phenomenon only lasts for a few pulse periods 



and not on a timescales of tens of days as detected for intermittent pulsars. A second argument that was proposed is that the intermittent 
phenomenology could be due to precession, which is an effect by which the pulsar undergoes a slow periodic wobble, thus moving the beams 
of radio radiation out of our line of sight. What ruled out this idea, however, is that precession certainly cannot produce a transition from the 
ON to the OFF state in less than 10s. 

A more convincing explanation proposed by Lvne 1 20091) and Gurevich & Istomir] 1 2007 ) is that there is a global failure of charge 
particle currents in the magnetosphere. In particular, the changes in the radio emission would be due to the presence or absence of a plasma 
whose current flow provides the expected extra torque on the star. In this model, the open field lines above the magnetic pole become depleted 
of charged radiating particles during the OFF states and the rotational slow-down, z>off, is produced by a torque dominated by magnetic 
dipole radiation. On the contrary, when the pulsar is ON, z>on is enhanced by an additional torque provided by the outflowing plasma. In 
other words, during the ON state the energy release is due to the current losses only, while during the OFF state it is due to the magnetodipole 
vacuum radiation (in this case, it is not the plasma-filled magnetosphere). In spite of its plausibility, this idea suffers from the lack of a 
physical mechanism for changing the plasma flow in the magnetosphere in such a drastic way. 

What we propose here is indeed an alternative idea based on our results about oscillating magnetospheres. As we have shown in Sec. [4] 
the energy loss by the pulsar can be significantly altered by the stellar oscillations. Therefore, it is reasonable to assume that during the 
ON state the stellar oscillations create relativistic wind of charged particles by virtue of the additional accelerating electric field. In a period 
of about lOd stellar oscillations are damped and the OFF period starts. The quasi-periodic glitch, whose driving mechanism is still largely 
obscure, is the only plausible excitation mechanism of oscillations for isolated pulsars and it would also be responsible for the emergence of 
new ON states with a certain periodicity. It is well known that the radio emission of pulsar is negligible in the overall energy budget, usually 
constituting a very tiny fraction of the pulsar spin-down rate, less then 10~ 3 . Most of the energy flux is carried away by the relativistic pulsar 
wind and does not reveal itself during pulsar emission. Therefore, switching radio emission on or off cannot change the pulsar spin-down 
significantly enough to be detected by observations. As a result, we propose that !>off ~ von and that pulsar quasiperiodic glitch is the real 
reason for i) the sudden increase of rotational energy at the end of the OFF state and for ii) the excitation of the stellar oscillations which 
switches pulsar radio e mission on. 

As suggested by Zhang etail 1 2007 ). the transition from the OFF to the ON state of intermittent pulsar would correspond to the 
reactivation of a dead pulsar above the 'death line' in the P — B diagram, thus becoming occasionally active only when the conditions for 
pair production and coherent emission are satisfied. It is wo rth stressing that it is very diffic ult to define an exact line in the P — B diagram 
for rotating neutron stars. However, in a recent investigation Ahmedov & Morozoval 1 2009h showed that oscillating but non rotating neutron 
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star remain below the death line for the majority of known radio pulsars. That result, when combined to the findings of the present study, 
suggest that oscillations in a rotating star could be the key ingredient for explaining the transient re-activation of intermittent pulsars. 

The connection that we have proposed here between the intermittent pulsar phenomenology and the presence of oscillations in the 
magnetosphere of rotating magnetized neutron stars will be exhaustively investigated in a more quantitative way in a forthcoming paper. 



6 CONCLUSIONS 

In this paper we have studied the astrophysical pr ocesses in the polar cap region of the magnetosphere of an oscillating neutron star. The 



background spacetime is given by the metric of lHartle & Thornd |l968) within the slow rotation approximation. The novelties of our 



analysis consists in quantifying the contributions of stellar oscillations in a general relativistic framework. In particular, we have computed 
the general-relativistic corrections to the Goldreich-Julian charge density, to the electrostatic scalar potential and to the component of the 
el ectric field paralle l to the magnetic field lines in the polar cap region when toroidal stellar oscillations are present. As already remarked 
by Timokhir] i 2007 ). the effective electric charge density i.e., the difference between the Goldreich-Julian charge density pgj (proportional 



to H • B for rotating stars and to u ■ B for oscillating stars in the flat space-time case) and the electric charge density (proportional to B) 
in the oscillating star magnetosphere is responsible for the generation of an electric field parallel to the magnetic field lines. Such difference 
vanishes only at the surface of the star while in general it becomes significant l y larg e at some distance r from the surface, due to the fact that 



p can not compensate pa. As already pointed out by Muslimov & Tsyganl 1 1992 ). general relativistic terms arising from the dragging of 
inertial frames give very important additional contribution to this difference. These terms depend on the radial distance from the star as 1/r 3 
and have important influence on the value of accelerating electric field generated in the magnetosphere near the surface of the neutron star. 

Our solutions for the accelerating electric field for the oscillating and rotating magnetized neutron stars may have some significant 
implications for pulsar polar cap models. These models assume that charged particles are accelerated above the polar caps, initiating pair 
cascades through one-photon pair creation of photons. The electric field induced by the stellar oscillations becomes therefore very important. 
Thus, the potential drop at the pair formation front, and the total energy gained by particles in the open field region is larger for the oscillating 
star. Since the contribution from the stellar oscillations to the electric field depends on the amplitude of stellar oscillations, pulsars having 
larger K = fj(l)/flR will have larger accelerating potential drops. Our main conclusions about oscillating and rotating neutron stars can be 
summarized as follows: 

(i) The oscillation regime of particle ejection from the stellar surface increases the total power carried away by relativistic primary particles 
relative to the purely rotating regime. Moreover, the fluctuation of the charge density of particles ejected from the stellar surface modulates 
the particle energy along a field line. 

(ii) The energy losses along the open magnetic field lines in the polar cap region and due to toroidal oscillations are significantly larger 
than the rotational energy losses for the m = 1 modes of oscillation. In particular, the energy losses of the mode (I, m) — (2, 1) can be a 
factor 8 larger than the rotational energy losses, even for an oscillation amplitude at the star surface as small as r\ = 0.05 Q R. 

(iii) The oscillation-induced inhomogeneity of the physical conditions at the stellar surface may substantially affect the global electrody- 
namics within the inner magnetosphere of a neutron star. 

(iv) The new dependence obtained for the energy losses on the oscillating behavior reflects in a new relation, namely Eq. d58l l. between 
the product PP and the amplitude of the oscillation at the star surface. In cases when the moment of inertia of the star is known with good 
accuracy, such a relation will allow to fully appreciate the effects of oscillations on pulsar magnetospheres. 

Finally, we have proposed a connection between the phenomenology of intermittent pulsars, characterized by the periodic transition 
from active to dead periods of radio emission in few observed sources, with the presence of an oscillating magnetosphere. In particular, 
we propose that, during the active state, star oscillations induced by periodic glitches of the neutron star create relativistic wind of charged 
particles by virtue of the additional accelerating electric field. After a timescale of the order of tens of days stellar oscillations are damped, 
and the pulsar shifts below the death line in the P — B diagram, thus entering the OFF invisible state of intermittent pulsars. This seminal 
idea, proposed here on a qualitative level, will be further explored in a future work. 
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